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Abstract 

We present X-ray emission characteristics of the massive 0-type stars DH 
Cep and HD 97434 using archival XMM-Newton observations. There is no 
convincing evidence for short term variabihty in the X-ray intensity during 
the observations. However, the analysis of their spectra reveals X-ray struc- 
ture being consistent with two-temperature plasma model. The hydrogen 
column densities derived from X-ray spectra of DH Cep and HD 97434 are 
in agreement with the reddening measurements for their corresponding host 
clusters NGC 7380 and Trumpler 18, indicating that the absorption by stellar 
wind is negligible. The X-ray emission from these hot stars is interpreted in 
terms of the standard instability-driven wind shock model. 

Keywords: star:X-ray - star:binary - star:wind - star:individual (DH Cep 
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1. Introduction 

Massive, hot and luminous 0-type stars have powerful stellar winds, with 
mass loss rates of the order of 10"^ to 10~^ M© yr~^ and velocities of several 
thousand km s~^. These winds play a crucial role both in the stellar evolution 
as well as the galactic evolution. The existence of strong winds is supposed to 
be a primary reason for X-ray emission in the early type stars. Study of the 
X-ray emission from 0-type stars is important for understanding the high 
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energy processes that generate hot gas in the wind or in the atmosphere of hot 
stars. X-ray emission from 0-type stars was first discovered scrcndipitously 
with the EINSTEIN observatory (Harndcn ct al. 1979). Most 0-type stars 
are reasonably bright (10^^ < Lx < 10^"^ erg s~^ ) and have a soft thermal 
spectrum with kT < 1 keV. The X-ray luminosity of these sources is found 
to scale with their bolometric luminosity as Lx/Lboi ~ 10~^ to 10~^ (Long 
& White 1980, Pallavicini et al. 1981, Singh & Naranan 1982, Bcrghofcr 
et al. 1997; Sana et al. 2006). It is generally thought to be produced 
as a result of shocks, with velocity jumps up to a few hundred km s~^ , 
generated throughout the stellar wind because of dynamic instabilities (Lucy 
& White 1980; Owocki & Cohen 1999; Kudritzki & Puis 2000). This scenario 
is now called a standard model for X-ray emission from early-type stars 
(Rauw 2008). Although, most of the X-ray emission is dominated by soft X- 
rays as expected from wind-shock model (e.g. Feldmeier, Puis & Pauldrach 
1997), a significant and an unexpected harder component is also seen in the 
X-ray spectra (e.g. Leutenegger & Khan 2003; Bross et al. 2007; Stelzer et 
al. 2005). High quality observational data from XMM- Newton and Chandra 
observatories are indeed crucial to provide a strong test of the latest models. 
In order to understand the X-ray emission processes from 0-type stars, we 
have analyzed XMM- Newton data for two 0-type stars DH Cep and HD 
97434. The results of our analysis are presented here. 

DH Cep is an 06 V -|- 07 V spectroscopic binary system and a member 
of the cluster NGC 7380 (distance = 3.7 kpc, age = 10 Myr; Underbill 
1969; Massey, Johnson & De Gioio-Eastwood 1995; Hilditch, Harries, & Bell 
1996). It has an orbital period of 2.11 d (Penny, Gies & Bagnuolo 1997). 
A low signal, Einstein X-ray observation is reported by Chlebowski, Harden 
& Sciortino (1989) with Lx ~ 5.37 x 10^^ erg s"^ in the energy band 0.2-3.5 
keV. 

HD 97434 is an 09 V star (Skiff 2007), with a membership probability 
of 78% (Dias et al. 2006) in the cluster Trumpler 18 (distance = 1.5 kpc; 
age = 15 Myr; Vazquez & Feinstein 1990). The Einstein X-ray data for 
the source are cataloged by Chlebowski, Harden & Sciortino (1989) with 
Lx < 3.4 X 10^^ erg s~^ (corrected for the distance of 1.5 kpc) in the energy 
band 0.2-3.5 keV. 
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2. XMM-Newton Observations and Data Reduction 

We present analysis of the archival data obtained with the XMM-Newton 
observatory. Table [1] summarizes the XMM-Newton observations. Our anal- 
ysis is based on CCD images, lightcurves and spectra from the European 
Photon Imaging Cameras (EPIC). Data were acquired simultaneously with 
EPIC-PN camera (Striider et al. 2001) and two nearly identical EPIC-MOS 
(MOSl and M0S2;Turner et al. 2001) cameras. Data reduction followed 
standard procedures using XMM-Newton Science Analysis System software 
(SAS version 8.0.0) with updated calibration files (Ehle et al. 2004). A de- 
tailed description of the data reduction and analysis procedure is given in 
Bhatt et al. (2010). Owing to poor counts below 0.3 keV and above 7.5 
keV, our analysis is limited to the energy band 0.3-7.5 keV. For the star 
HD 97434, the observations are not made with the MOSl detector, therefore 
we used only the data from PN and M0S2 detectors for further analysis. 
For the star DH Cep, the observations have been taken during the binary 
phase of 0.08-0.25. The binary phases of DH Cep were determined using the 
ephemeris HJD = 2441905.805+2. 110932E (Hilditch, Harries, & Bell 1996; 
Penny, Gies & Bagnuolo 1997). 

Source photons were extracted from circular region of 20" around the 
source position to generate the light curves and spectra of the source DH 
Cep and HD 97434. Background was estimated from a number of empty 
regions close to the X-ray source in the detector. X-ray spectra of the sources 
were generated using SAS task especget, which also computes the photon 
redistribution as well as the ancillary matrix. Finally, the spectra were re- 
binned to have at least 20 counts per spectral bin for both the sources. 

3. Analysis 

Figure [1] shows the background and the background subtracted X-ray 
lightcurves of the stars DH Cep and HD 97434. We have binned the event 
lists into 1000 and 1500s time intervals for the stars DH Cep and HD 97434, 
respectively. We have performed test to measure the significance of the 
deviations from the mean count rates. The probability of the variability 
(Pvar) to quantify the constancy of the data over the time scale of observations 
was calculated. The results are displayed in the respective panels of Figured] 
The value of Pvar was calculated to be > 99% for the star DH Cep. However, 
it was calculated to be < 50% for the star HD 97434. This indicates that 
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the light curve of the star DH Cep shows significant variabihty. We further 
analysed the light curves of DH Cep and found that the observed variability 
may arise due to a few outlier points. If we neglect the four outlier points 
of the light curve of DH Cep, the probability of variability decreases to 80%. 
Therefore, we have plotted the 3a limit in Figure [1] represented by a shaded 
region around the constant count rate, i.e., mean value of count rate for the 
data sample. Here, a represents the standard deviation in count rates. This 
analysis shows that none of the data points in any of the light curves were 
lying outside of 3a limit. Further, the variability trend seen in the PN data 
was not correlated with the MOS data and all the data points in the light 
curves were within 3a limits. Therefore, it appears that the significance of 
the variability in the light curves is very poor and cannot be detected by the 
present data set. 

X-ray spectra of DH Cep and HD 97434 are shown in Figure [2l Emis- 
sion lines of Fe XVII (0.8 keV), Ne X (1.02 keV), Mg XII (1.47) and Si 
XIII (1.853 keV) are seen in their respective spectra. X-ray emission from 
massive stars is thought to be produced from thermal emission, generated 
from shock-heated plasma. These hydrodynamic shocks may occur in the 
stellar winds from massive stars. However, in massive binaries they may oc- 
cur in wind- wind collision zone (Bhatt et al. 2010 and references therein). 
Therefore, we have modeled the spectra using (a) plane-parallel shock model 
(PSHOCK; Borkowski, Lyerly & Reynolds 2001), and (b) models of Astro- 
physical Plasma Emission Code (apec; Smith et al. 2001), as implemented 
in the XSPEC version 12.3.0. A ^ minimization gave the best fitted model 
to the data. We corrected for the local absorption in the line-of-sight to the 
source (Nh) using the photoelectric absorption cross sections according to 
Balucihska-Church & McCammon (1992) and modeled it as phabs (photo- 
electric absorption screens). The abundances in the absorber were allowed to 
vary during the fitting procedure The abundances were varied with respect 
to solar abundances adopted from Lodders (2003). 

First, we fitted PSHOCK plasma model to derive their spectral features. 
The constant temperature PSHOCK plasma model was considered without 
incorporating mass-loss and orbital parameters of massive stars. However, 
the model does account for non-equilibrium ionization effects and assumes an 
equal electron and ion temperature. The best fit PSHOCK models to the data 
are shown in left and right top panels of Figure [2] and the best fit parameters 
are given in Table |2j The best-fit X-ray temperatures deduced from X-ray 
spectra are found to be 0.64+°;°^ keV and 0.64^°;°^ keV for DH Cep and HD 
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97434, respectively. 

Second, the simplest spectral APEC model was used for the spectral fit- 
ting. Here, we first considered the case of an isothermal hot ionized gas as 
the one temperature (IT) APEC model. This model is expressed as phabs 
X APEC. The IT APEC models with either solar or subsolar photospheric 
abundances were found to be unacceptable for both the stars due to high 
value of xl- Therefore, for the next level of complexity, we used models with 
two temperatures (2T) expressed as phabs (apec+apec). In 2T models the 
first component represents the "cooler" ionized gas and the second compo- 
nent represents the "hotter" ionized gas. The 2T models were found to be 
acceptable but required subsolar abundances, i.e., Z ^ 0.2 Zq. The subsolar 
abundances have also been reported for 15 other massive stars by Zhekov 
& Palla (2007). The best fit 2T models are displayed as histograms along 
with the data in the left and the right bottom panels of Figure |2] for the 
stars DH Cep and HD 97434, respectively. The best-fit parameters for these 
models were obtained by x^-minimization technique and are given in Table [2] 
with their corresponding xl- We noted that 2T APEC model fitted these 
spectra better than pshock model with a smaller value of xl (see Figure [2] 
and Table 12]). The cooler and hotter X-ray temperatures deduced from X-ray 
spectra are found to be respectively 0.62^q;q2 keV and > 1.89 keV for DH Cep 
and 0.32+°:?^ keV and 0.74+°;^^ keV for HD 97434, respectively. The best fit 
values of Nh for DH Cep and HD 97434 are in agreement with the values 
of Nh derived from optical data for NGC 7380 (0.32 x 10^2 cm'^ ; Massey, 
Johnson & De Gioio-Eastwood 1995), and Trumpler 18 (0.14 x 10^^ cm"^ 
; Vazquez & Feinstein 1990), respectively. The Nh from optical data was 
estimated using the relation, Nh = 5.0 xlO^^ x E(B — V) cm"^ (Vuong et al. 
2003), where E(B — V) = Av/3.1, assuming a normal interstellar reddening 
law towards the direction of the cluster. The extra absorbing hydrogen col- 
umn density in addition to the interstellar contribution was not needed in 
the spectral fits implying that the absorption by the stellar wind is negligible 
(see also Zhekov & Palla 2007; Sana et al. 2006). 

4. Discussion 

4.I. X-ray Variability 

The wind-shock model predicts a significant short-term variability in the 
X-ray fiux as the shocks fade and grow on short time scales. However, it 
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has not been observed so far. We have also explored the possibility of short- 
term X-ray variability of DH Cep and HD 97434 over the observation time 
span (<40 ks). The timing analysis reveals a lack of short term variability 
in the X-ray intensity for both the massive stars. The lack of short term 
variability could be due to the wind fragmentation, so that individual X-ray 
fluctuations are smoothed out over the whole emitting volume, leading to a 
rather constant output ( Feldmeier, Puis & Pauldrach 1997). Several other 
massive stars, e.g., HD 159176 (De Becker et al. 2004) and HD 47129 (Linder 
et al. 2006) have also shown lack of variability in their light curves. 

The longer-term variability can not be probed using present data. The 
long-term variability (few days) from 0-type stars was reported by Snow, 
Cash & Grady (1981) for the first time, but their results were based on the 
low sensitivity detectors available at that time. No significant variability was 
detected by the ROSAT All- Sky Survey in the X-ray emission of 57 OB stars 
over a time-scale of ~2 d (Bcrghofer & Schmitt 1995). 

4-2. X-ray temperatures of plasm,a 

Our analysis of XMM-Newton data reveals reveal two-component X-ray 
temperature structure for massive stars. The temperatures corresponding to 
the cool and the hot components of HD 97434, and the cool component of DH 
Cep are found to be less than 1 keV. However, the temperature corresponding 
to the hot component of DH Cep was found to be more than 1.89 keV. 
The derived values of temperatures are similar to the 0-type stars in the 
NGC 6231 cluster (Sana et al. 2006) and in the Carina OBI association 
(Antokhin et al. 2008). Generally for 0-type stars, the best fits to good 
quality data can usually be achieved by the sum of two thermal components 
at about 0.3 and 0.7-1 keV (Giidel & Naze 2009). The observed soft X-ray 
emission could originate from radiation-driven instabilities in stellar winds 
(Lucy 1982). The wind-shock model predicts the intrinsic instability of the 
line driving mechanism. Indeed, the velocity in an unstable wind is not the 
same everywhere and the fast-moving parcels of material will overcome the 
slow moving material, generating shocks between them. This model estimates 
the shock velocities (vghock) by the relation (Lucy 1982, Luo, McCray & 
Maclow 1990) 




(1) 
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Adopting the mean particle weight /i ^^^0.62 for 0-type stars (Cassinelh et 
ah 2008) we derive the "average" value of shock velocities correspond to the 
cool components (kTi) 716l]^2 and 514]|^^q2 km s~^ for the stars DH 

Cep and HD 97434, respectively. However, the shock velocities correspond to 
the hot components (kT2) are found to be > 1563 km s'^ and 782+^^ km s'^ 
for DH Cep and HD 97434, respectively. These values are about a factor of 2 
larger than the predicted by radiative shock model of Lucy (1982). However, 
the evolved version of the standard model by Owocki, Castor & Rybicki 
(1988) predicts X-ray emission upto 1 keV from wind shocks model. The 
temperature corresponding to the hot component of DH Cep is found to be 
more than 1.89 keV. It may arise from wind collision zone of the binary 
system (see Bhatt et al. 2010 and references therein), but we are not able 
to draw any firm conclusion on the basis of it alone. Therefore, it appears 
that the cool as well as hot temperature components from DH Cep and HD 
97434 could be generated by instabilities in radiation-driven wind shocks. 

4-3. X-ray Luminosity 

For massive stars DH Cep and HD 97434, the ratios of X-ray to bolometric 
luminosities, log(Lx/Lboi ), are found to be -6.7 and -7.3 (see Table [2] ) in 
the energy band 0.3-7.5 keV, respectively, which are broadly consistent with 
the relation derived for similar kind of 0-type stars (Sana et al. 2006). 

X-ray luminosities derived from XMM-Newton data are found to be an 
order of magnitude lower than the Lx derived from Einstein observations for 
both the massive stars (see Table [2] and 51]) • These discrepancies can not 
be explained by difference in sensitivity ranges of the various instruments. 
A similar results have been found by De Becker et al. (2004) for the star 
HD 159176 after the comparison of XMM-Newton, Einstein and ROSAT 
data. We found the X-ray luminosity 4.82 x 10^^ and 1.95 x 10'^^ in energy 
range 0.3-2.0 keV for the DH Cep and HD 97343, respectively, using XMM- 
Newton data. The choice of the energy range is nearly a similar to the energy 
range in Einstein observations, i.e., 0.2-3.5 keV. We further converted the 
Einstein IPC count rates into flux using the WebPIMMdll, where we assumed 
a thermal plasma model with temperature of 0.5 keV, and the hydrogen 
column density of 3.9 x lO^^ for DH Cep and 1.8 x lO^^ for HD 97434 (see 
Table 2). The Einstein IPC count rates are 0.0166 and < 0.0103 counts s~^ 



^http://heasarc. gsfc.nasa.gov/Tools/w3pinims. html 
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for the stars DH Cep and HD 97343, respectively (Chlebowski et al. 1989). 
The X-ray himinosity thus calulated was 2.3 x 10^^ and < 1.2 x 10'^^ for the 
stars DH Cep and HD 97434, respectively. These value are consistent with 
that obtained from XMM-Newton obsertions. Therefore, it appears that the 
model used by Chlebowski et al. (1989) to convert the Einstein count rates 
into luminosities may be responsible for these discrepancies. 

5. Summary and Conclusions 

Using the archival XMM-Newton observations, we present, for the first 
time, the temporal and spectral analysis of X-ray emission from 0-type stars 
DH Cep and HD 97434. The main results are as follows: 

1. There is no firm evidence for short term variability in the X-ray inten- 
sity during the observations time span (< 40 ks) for both the stars. 

2. The values of Nh derived from the best fit X-ray spectra are consistent 
with the optical estimates of reddening. 

3. X-ray spectra of these stars are fitted consistently with two-temperature 
plasma models. The X-ray emitting plasma is found to be generated at 
a temperature of lesser than 1 keV. It could originate from small shocks 
in the radiation-driven outflows. A hotter component is indicated in 
DH Cep. 

4. The best-fit values of abundances are found to be 0.2 times of the solar 
abundances. 
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Table 1: Journal of XMM-Newton Observations of the objects. 



Object Name 



DH Cep 



HD 97434 



Observation ID 
Exposure Time (s) 
Start time UT (hh:mm:ss) 
Usable time(ks) (MOSl)t 
Usable time(ks) (M0S2) 
Usable time(ks) (PN) 
EPIC filter 
Offset from Target (arcmin) 



0205650101 
31413 
19 Dec 2003 02:02:12 
25.23 
25.46 
19.52 
Thick 
0.000 



0051550101 
40822 
06 Feb 2002 01:13:19 

39.20 
35.85 
Medium 
3.68 



t : The observations have not done in Prime Full window mode for MOSl 
detector. 



Tabic 2: The best-fit spectral parameters of the stars DH Cep and HD 97434. 



Parameters 




DH Cep 


HD 97434 






PSHOCK 


2T APEC 


PSHOCK 


2T APEC 


Nh(1022 cm-2) 




42+°-°^ 


n Qo+u.ua 


-0.04 


n 18+0.05 


kTi(keV) 






62+°-°2 


'-'■"^-0.04 


32+°-°^ 

'-'•'-'^-0.17 


kT2(keV) 






> 1.89 




74+0-13 

'^-0.12 


Z/Ze 






0.22^°f 




21+°-^^ 


t(10^^ s cm-3) 




15.14™ 




9mtl',f 




EMi(105^ cm-3) 




240.26li:«^ 








EMailO^^ cm-3) 










o /1O+2.02 


X^/d.o.f. 




1.32(181) 


1.18(180) 


1.36(75) 


1.31(74) 


Fx(10~^^ erg cm" 


-2 s-i) 


0.33 


0.34 


0.08 


0.08 


log Lx(erg s"^) 




32.73 


32.75 


31.31 


31.32 


log(Lx/Lboi)[0.3- 


-7.5] 


-6.73 


-6.71 


-7.29 


-7.29 


log(Lx/Lboi)[0.3- 


-2.0] 


-6.77 


-6.77 


-7.31 


-7.31 


log(Lx/Lboi)[2.0- 


-7.5] 


-7.76 


-7.55 


-8.56 


-8.48 
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Figure 1: X-ray light curves for the stars DH Cep (left panel) and HD 97434 (right panel). 
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Figure 2: Best-fit X-ray spectra of massive stars, left panel: for DH Cep. right panel : 
for HD 97434. The solid line histograms represent the best fit pshock (in upper panels) 
and 2T APEC model (in lower panels) with absorption model PHABS for PN, MOSl and 
M0S2 data. 
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